A new Monte Carlo muon generator for cosmic-ray muon applications

Cosmic rays, thanks to their ubiquity and high penetration capability, have been successfully used in scientific research ever since their discovery. As soon as their knowledge improved, applications in the civil/environmental field were also developed: muon radiography (or muography, based on the flux attenuation) and muon tomography (based on the scattering angle) have been used to study the inner structure of volcanoes, to seek hidden rooms in Egyptian pyramids, to search for heavy metals in containers, and so on. And besides these imaging techniques, cosmic ray muons are also widely used for detector testing and alignment practically in every Nuclear Physics or Particle Physics experiment. Since most of these applications are sensitive to the angular and momentum distribution of cosmic muons, an accurate modelling of these distributions is a key feature for any generation tool conceived to simulate the cosmic muon flux. This can make the generator quite time-consuming, which is a strong limit when one needs to reach high statistics or to study large structures. A new Monte Carlo generator for cosmic-ray muons, named Efficient COsmic MUon Generator (EcoMug for short), especially designed to be fast (≳ 10^5 muons generated per second on a standard machine) without losing accuracy, is presented here. It is written as a header-only C++11 library, ready to be integrated into whatever C++ code, in particular C++ code based on Geant4 simulation tool. By default, EcoMug relies on a simple and effective parametrisation of the experimental data of cosmic ray differential flux at sea level, taken from the literature, but the library is written in such a way that every user can easily replace it with his own user-defined parametrisation. Unlike other tools, EcoMug is able to generate muons from different kind of surfaces (plane, cylinder and half-sphere), while keeping the correct angular and momentum distribution of generated tracks inside a fiducial volume. This allows to optimise the generation surface according to the system under study, and leads to a further improvement of the overall simulation efficiency. In this contribution we will present the main features of EcoMug, starting from its mathematical foundation, and eventually showing some interesting applications

Noise reduction in muon tomography for detecting high density objects

The muon tomography technique, based on multiple Coulomb scattering of cosmic ray muons, has been proposed as a tool to detect the presence of high density objects inside closed volumes. In this paper a new and innovative method is presented to handle the density fluctuations (noise) of reconstructed images, a well known problem of this technique. The effectiveness of our method is evaluated using experimental data obtained with a muon tomography prototype located at the Legnaro National Laboratories (LNL) of the Istituto Nazionale di Fisica Nucleare (INFN). The results reported in this paper, obtained with real cosmic ray data, show that with appropriate image filtering and muon momentum classification, the muon tomography technique can detect high density materials, such as lead, albeit surrounded by light or medium density material, in short times. A comparison with algorithms published in literature is also presented

Precision measurements of Linear Scattering Density using Muon Tomography

We demonstrate that muon tomography can be used to precisely measure the properties of various materials. The materials which have been considered have been extracted from an experimental blast furnace, including carbon (coke) and iron oxides, for which measurements of the linear scattering density relative to the mass density have been performed with an absolute precision of 10%. We report the procedures that are used in order to obtain such precision, and a discussion is presented to address the expected performance of the technique when applied to heavier materials. The results we obtain do not depend on the specific type of material considered and therefore they can be extended to any application.Comment: 16 pages, 4 figure

Search for charginos in e+e- interactions at sqrt(s) = 189 GeV

An update of the searches for charginos and gravitinos is presented, based on a data sample corresponding to the 158 pb^{-1} recorded by the DELPHI detector in 1998, at a centre-of-mass energy of 189 GeV. No evidence for a signal was found. The lower mass limits are 4-5 GeV/c^2 higher than those obtained at a centre-of-mass energy of 183 GeV. The (\mu,M_2) MSSM domain excluded by combining the chargino searches with neutralino searches at the Z resonance implies a limit on the mass of the lightest neutralino which, for a heavy sneutrino, is constrained to be above 31.0 GeV/c^2 for tan(beta) \geq 1.Comment: 22 pages, 8 figure

Search for composite and exotic fermions at LEP 2

A search for unstable heavy fermions with the DELPHI detector at LEP is reported. Sequential and non-canonical leptons, as well as excited leptons and quarks, are considered. The data analysed correspond to an integrated luminosity of about 48 pb^{-1} at an e^+e^- centre-of-mass energy of 183 GeV and about 20 pb^{-1} equally shared between the centre-of-mass energies of 172 GeV and 161 GeV. The search for pair-produced new leptons establishes 95% confidence level mass limits in the region between 70 GeV/c^2 and 90 GeV/c^2, depending on the channel. The search for singly produced excited leptons and quarks establishes upper limits on the ratio of the coupling of the excited fermio

Hadronization properties of b quarks compared to light quarks in e+e- -> q qbar from 183 to 200 GeV

The DELPHI detector at LEP has collected 54 pb^{-1} of data at a centre-of-mass energy around 183 GeV during 1997, 158 pb^{-1} around 189 GeV during 1998, and 187 pb^{-1} between 192 and 200 GeV during 1999. These data were used to measure the average charged particle multiplicity in e+e- -> b bbar events, _{bb}, and the difference delta_{bl} between _{bb} and the multiplicity, _{ll}, in generic light quark (u,d,s) events: delta_{bl}(183 GeV) = 4.55 +/- 1.31 (stat) +/- 0.73 (syst) delta_{bl}(189 GeV) = 4.43 +/- 0.85 (stat) +/- 0.61 (syst) delta_{bl}(200 GeV) = 3.39 +/- 0.89 (stat) +/- 1.01 (syst). This result is consistent with QCD predictions, while it is inconsistent with calculations assuming that the multiplicity accompanying the decay of a heavy quark is independent of the mass of the quark itself.Comment: 13 pages, 2 figure

Search for supersymmetric particles in scenarios with a gravitino LSP and stau NLSP

Sleptons, neutralinos and charginos were searched for in the context of scenarios where the lightest supersymmetric particle is the gravitino. It was assumed that the stau is the next-to-lightest supersymmetric particle. Data collected with the DELPHI detector at a centre-of-mass energy near 189 GeV were analysed combining the methods developed in previous searches at lower energies. No evidence for the production of these supersymmetric particles was found. Hence, limits were derived at 95% confidence level.Comment: 31 pages, 14 figure

Search for lightest neutralino and stau pair production in light gravitino scenarios with stau NLSP

Promptly decaying lightest neutralinos and long-lived staus are searched for in the context of light gravitino scenarios. It is assumed that the stau is the next to lightest supersymmetric particle (NLSP) and that the lightest neutralino is the next to NLSP (NNLSP). Data collected with the Delphi detector at centre-of-mass energies from 161 to 183 \GeV are analysed. No evidence of the production of these particles is found. Hence, lower mass limits for both kinds of particles are set at 95% C.L.. The mass of gaugino-like neutralinos is found to be greater than 71.5 GeV/c^2. In the search for long-lived stau, masses less than 70.0 to 77.5 \GeVcc are excluded for gravitino masses from 10 to 150 \eVcc . Combining this search with the searches for stable heavy leptons and Minimal Supersymmetric Standard Model staus a lower limit of 68.5 \GeVcc may be set for the stau mas

Measurement of the t t-bar production cross section in the dilepton channel in pp collisions at sqrt(s) = 7 TeV

The t t-bar production cross section (sigma[t t-bar]) is measured in proton-proton collisions at sqrt(s) = 7 TeV in data collected by the CMS experiment, corresponding to an integrated luminosity of 2.3 inverse femtobarns. The measurement is performed in events with two leptons (electrons or muons) in the final state, at least two jets identified as jets originating from b quarks, and the presence of an imbalance in transverse momentum. The measured value of sigma[t t-bar] for a top-quark mass of 172.5 GeV is 161.9 +/- 2.5 (stat.) +5.1/-5.0 (syst.) +/- 3.6(lumi.) pb, consistent with the prediction of the standard model.Comment: Replaced with published version. Included journal reference and DO